arXiv: 1503.03953v4 [quant-ph] 8 Oct 2015 


Large-N scaling behavior of the quantum fisher information in the Dicke model 
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Quantum Fisher information (QFI) of the reduced two-atom state is employed to capture the 
quantum criticality of the superradiant phase transition in the Dicke model in the infinite size and 
finite-systems respectively. The analytical expression of the QFI of its ground state is evaluated 
explicitly. And finite-size scaling analysis is performed with the large accessible system size due to 
the effective bosonic coherent-state technique. We also investigate the large-size scaling behavior of 
the scaled QFI of the reduced A^-atom state and show the accurate exponent. 
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I. INTRODUCTION 


Quantum Fisher information (QFI), one of the quan¬ 
tum information-based tools, is a basic concept in quan¬ 
tum estimation theory, which depicts the theoretical 
bound for the variance of an estimator [IH3. The Fisher 
information is the central notion in parameter estima¬ 
tion due to the Cramer-Rao inequality, which sets a ba¬ 
sic lower bound to the variance of any unbiased estima¬ 
tor in terms of the Fisher information Q. Consider the 
problem of estimating a unknown parameter 9 from a 
quantum state p{9). The value of 9 can be estimated 
from the measurement results of a proper physical ob¬ 
servable. Since the precision of the estimation is limited 
by unavoidable measurement errors, the inverse of the 
QFI provides the lower bound of the error of the estima¬ 
tion of 9. 

Latterly, a new emphasis has emerged in which QFI 
is related to properties of interacting many-body sys¬ 
tems. This approach is being pursued most vigorous^ 
in connection with quantum phase transition (QPT) [9j, 
as it is hoped that the QFI may shed light upon the 
dramatic effects occurring in critical systems. QPT and 
quantum-critical phenomena occur at zero temperature 
in many-body quantum systems. A dramatic change of 
order parameters exhibits at the critical point, which are 
induced by the change of parameters in quantum criti¬ 
cal systems. There is much on-going interest to test the 
quantum criticality of the QFI in the proximity of phase 
transitions (l0l - [l^ . It is expected to characterize the sin¬ 
gularity of the QPT from quantum estimation perspec¬ 
tive by driving the system toward critical points. In this 
framework, we consider the scaling exponents of the QFI 
at the critical point in the Dicke model 0 , which is a 
well-known quantum collective atoms model. 

The Dicke model describes the interaction of N two- 
level atoms with a single bosonic mode. The QPT was 
explored in the Dicke model, exhibiting a superradiant 
phase transition in the thermodynamics limit 
Although the Dicke model cannot be solved analytically, 
an extended bosonic coherent state approach can solve 
the Dicke model numerical accurately for large size sys¬ 
tems 0. For finite-size atoms Dicke model has been 
characterized in terms of entanglement of its ground 
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of QFI has not been well analyzed, except preliminarily 
results for the QFI of the field mode and A-atom state 
in the ground state [2^. A convincing finite-size scaling 
behavior of the QFI is still lacking. To the best of our 
knowledge, the finite-size studies are limited to numeri¬ 
cal diagonalization in the bosonic Fock state in small-size 
systems A < 35 [H, . Our paper is intended to pro¬ 

pose the QFI of the reduced two-atom state to study the 
quantum criticality of the QPT in the Dicke model and 
the finite-size scaling exponents. 

In this paper, we study the QFI of the reduced two- 
atom state in the Dicke mode in infinite size and finite-A 
systems respectively, giving the accurate finite-size scal¬ 
ing exponents. The paper is organized as follows. In Sec. 
II we review the definition of the QFI and its physical sig¬ 
natures by the parameter estimation theory. In Sec.Ill 
we introduce the background of the Dicke model, and de¬ 
fine the QFI of the reduced two-atom state to capture the 
quantum criticality in the large-A atoms system and in 
the thermodynamics limit. In Sec.IV the large A scaling 
behavior of both of the QFI of the two-atom state and A- 
atom state are calculated by the bosonic coherent-state 
technique, giving the accurate scaling exponents. Finally, 
we summarize our work in Sec.IV. 


II. GENERAL FORMALISM FOR THE QFI 

To begin with, we briefly review the parameter esti¬ 
mation theory and the QFI, which is applied to evaluate 
bounds of the variance of estimator for a parameter that 
we can get from a quantum state. To measure the preci¬ 
sion of the estimator 9, we consider a quantum state p(9). 
The generalized quantum Fisher information (QFI) F is 
defined as 


FQ{p{9))=rr[p{9)L\ (1) 

The symmetric logarithmic derivative operator L is de¬ 
termined by 


dep{9) = ]^[Lp{9) + p{9)L]. (2) 
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Assume that the spectral decomposition of the density 
operator is given by p{9) = with eigen¬ 

values Pi and eigenvectors \ipi) of p(0). And L can be 
solved by rewriting the above equation under the eigen- 
basis of piO). Then the QFI obtained can be written as 


Fq{.p{^)) = '^'^Pi(.{dePi\dePi) -\{p^\dePi)\'^) 

1^^ j|2 

The value of parameter 0 can be estimated through mea¬ 
suring p{9). From the QFI, we obtain the lower bound 
of the variance of the estimator for the parameter 9, 
given by the quantum Cramer-Rao (QCR) theorem: [1,0 
(69)'^ > 1/Fq{p{9)). From this inequality, it is obvious 
that the variance of the estimation is small for a large 
value of Fq{p{9)). 

The estimation of parameters in general quantum 
metrology process consists of the following steps. First, 
we prepare a quantum state Then the system un¬ 
dergoes the 0-dependent process with the phase- 

shift generator U, and evolves to the state p{9) = 
^-i0Up^^^ieu^ Finally, we estimate the parameter 9 
by measuring p{9). The variance of the estimation is 
bounded by the inverse of the QFI. As the eigenvalues of 
p{9) and pm are the same, the expression of QFI ( [31) is 
simplified as 


turn form BUS 


2A 

H = ujo^a + AJz H—+ a)J., 


where and a are the bosonic annihilation and creation 
operators of the single-mode cavity, A and w are the 
transition frequency of the two-level atom and the fre¬ 
quency of the single bosonic mode, A is the coupling 
constant. Jx and are the collective atomic oper¬ 
ators as Ja = = 3:, j/, z). And the 

Hilbert space of this algebra is spanned by the Dicke 
state {|j,m),m = -j,-j + 1,..., j - 1, j} with j = A/2, 
which is the eigenstate of and Jz- 

The Hilbert space of the total system can be expressed 
in terms of the basis {|n) ® |j,m)}, where |n) is the 
Fock state containing n bosons. To our best knowl¬ 
edge, it is very difficult to obtain convergent results for 
large number of atoms based on usual basis of the Fock 
states [11,120. An extended coherent state technique [13 
has been proposed to the Dicke model up to large N- 
atom system, which has been confirmed recently by com¬ 
paring with the results in terms of basis of the Fock 
states [11]. In the extended coherent state approach, 
the wave function can be expressed in terms of the basis 
{\Pm)}) ® |j,m)}, where is the bosonic extended 

coherent state 


Nt, 


k=0 


= E |0), , (5) 


Fq{p{ 9),U) =A^p,{SUf - ^ .^^\{p,\U\p,)f. 

■ 1 Pi ' Pj 

(4) 

where {SU)"^ = {(pi\U‘^\(pi) — \{(pi\U\pi)\‘^. For a pure state 
Pin = \p){'p\^ and the QFI is simplified as Fq{p{9), U) = 
^Y^l^iPi{5UY. It is obvious that the QFI depends on 
the quantum state pin and the choice of the phase-shift 
generator U. In our work, since the Dicke model is the 
collective atom model, the quantum state pin is prepared 
in a specific class. For collective models, all two-level 
atoms are completely equivalent and the ground states 
of collective models are invariant under the permutation 
group. We discuss the QFI of the reduced two-atom state 
which is a symmetric state and permutation-invariant 
in the Dicke model. The QFI is expected to capture 
the quantum criticality and signal the presence of QPT. 
Moreover, it is significant to compute the finite-size scal¬ 
ing exponent of the QFI for the universality of the QPT. 


III. THE QFI OF THE REDUCED TWO-ATOM 
STATE 

We study the QFI and its scaling behavior for N two- 
level atoms system in the Dicke model. The Dicke Hamil¬ 
tonian can be written in terms of the collective momen- 


where pm = 2m\/ (uj'/N), Ntr is the truncated bosonic 
number in the space of the new operator Am = a + Pm, 
|0)^ is the vacuum as ajO)^ = 0, and the coefficient Cm,k 
can be determined through the exact diagonalization. 
Then the ground state for the finite-A system takes the 
form of 

j 

1 '^)= \Prn)b(^\j,rn). ( 6 ) 

m=-j 

In the ground state it undergoes a transition from the 
normal to the super-radiant phase when increasing the 
coupling A through! a critical value of Ac = •\/wA/2 [H- 
[13, 0, IS, [30. In the super-radiant phase, the atomic 
ensemble spontaneously emits with an intensity propor¬ 
tional to A^ rather than A. 

Since there is on going interest in studying the connec¬ 
tion between the QFI and QPT. The QFI is expected to 
shed light on the dramatic effects occurring at the crit¬ 
ical point, providing to signal the presence of QPT. It 
is necessary to study the QFI in the ground state of the 
Dicke model. 

In the quantum metrology process, we choose the re¬ 
duced two-atom state p associated with arbitrary two 
atoms as a prepared state, which is obtained by trac¬ 
ing out the density matrix over all other A — 2 atoms 
and field mode. For collective atom ensemble, the re¬ 
duced two-atom state is particularly well suited since it 




3 


does not depend on the two atoms selected, all atoms 
being completely equivalent. Since the iV-atom state is 
the Dicke state \j, m >, which is a symmetric state. The 
reduced two-atom state can be extracted from the sym¬ 
metric multi-atom Dicke state \j,rn >. Due to the sym¬ 
metry of the state of A^-atoms under exchange of atoms, 
the reduced two-atom state p is invariant under the per¬ 
mutation Group. The matrix elements of the reduced 
two-atom state p can be expressed in terms of the aver¬ 
ages of the collective atomic operators, which has been 
addressed by Wang and J.Vidal in Ref. [sll, Then 
the reduced two-atom state can be given in the basis of 
two atoms {|U), lit), lU), Itt)} (with cr^|t) = |t) and 
<7z\i) = -|i)) as 


( V+ 

X*i_ X*J^ 

* 

u 

x+ 

w 

y 


x+ 

y 

w 

x_ 

V U 

X- X- 

V- 


(7) 


where the matrix elements can be represented by the ex¬ 
pected values of the collective spin operators 


v± 


x± 


w 


u 


-2N + 4(J2)±4(V - 

47V(iV- 1) ’ 

(jV-l)(J+)±([J+,J.]H-) 

2N{N-l) 

- 4(Jf) _ {Jl + Jl) - N/2 

47V(iV-l)’^ N{N-l) 

{Jl) 

N{N -\y 


where [A, B]+ = AB + BA, and w = y ior J2a=x,y,z = 

= t(t + 7). (Jc) and {Ja){a = x,y,z,+) mean 
the averages of the collective atomic operators over the 
ground state |G), Eq.( (B]). Thus, we can calculate the 
expected values of the collective atomic operators to de¬ 
termine the elements of the reduced two-atom state. 

Since there is a conserved parity If in the Dicke model, 
such that [i?, n] = 0, which is given by If = with 
the excitation number N = a''^a + Jz + N/2. The parity 
n possesses two eigenvalues ±1, depending on whether 
the number of quanta is even or odd. Then the Hilbert 
space of the total system is split into two noninteracting 
subspaces, resulting (G|J±|G) = 0. Then it is easily to 
find cci = 0 in Eq. ( B- Hence the reduced two-atom 
state can be shown in X form as 


/ v+ 0 0 u* 

_ 0 w y 0 

^ ~ 0 y w 0 

\u 0 0 

The reduced two-atom state facilitates the analytical 
evaluation of the QEI of the two-atom state in the fol¬ 
lowing. 

We consider an estimation of the parameter 9 in¬ 
troduced by the following unitary transformation S = 



exp(—i9a/) with the phase-shift generator erf = cr^ 0 
I 1^. Here I is the 2x2 identity matrix and az is the 
pauli matrix. From the definition in Eq. (|4|), the analyt¬ 
ical expression of the QEI are obtained as 




8piPj 

yy. Pi + Pj 




= 16( 




v+ +V- 2 


( 10 ) 


which are evaluated in detail in Appendix. We calculate 
the QEI in the infinite and finite size systems as Fq^oo 
and Fq^n respectively . 

We begin to discuss the QEI of the reduced two-atom 
state in X form by evaluating the expected values of ma¬ 
trix elements in the thermodynamics limit, in which the 
number of atoms becomes infinite. In this limit analysis, 
we first apply the Holstein-Primakoff transformation to 
change the collective ang ular oper ators to the boso n op¬ 
erators b{by by J+ = b^VN — b'^b, J = y/N — W bb, and 
= yb— N/2, where [&, 6'^] = 1 [l^. Then the displace¬ 
ments of the boson operators are introduced to depict the 
behaviors of super-radiation phase as -I- yfNa and 

= b^ — yfN/d. By means of the boson expansion ap¬ 
proach, we expand the Hamiltonian with respect to the 
new operator and as power series in 1 /N, 


H = NHo + N^^^Hi + ■■■, (11) 


where i?o = — dXa/dy/l — -f A(/3^ — 1/2) and 

Hi = -aoj + 2\I3{1 - /3V2)(ct -b c) -f A^(dt + d). By 
using large N expansions of H up to the 1/V, we obtain 
the ground state energy Ea{a,f3) as 

^£^^=u;a^-4XapVl^+A{/3^-^). (12) 

Minimizing the ground state energy gives 


oja-2XI3y/r^ = 0 (13) 

2aAyr^-^^-/3A = 0. 

Vt - P 

then we have 

= max{0, i(l - ^)}, (14) 

a = —/3y/l - /32, 

OJ 

where /r = 1 in the normal phase and p = (Ac/A)^ in the 
superradiant phase with the critical point Ac = \fijjK/2. 
Next we can derive the matrix elements of the reduced 
two-atom state p in Eq. ([9]) up to 0(1) 

V+ = V- = {1- (15) 

w = y = p\l-fi^), 

u = 





















4 


In the thermodynamics limit TV —>■ oo, the QFI of the 
reduced two-atom state Fq ^o can be expressed as 

^4 + (1_^2)2- (6) 

We know that I3^ + {1 — > 2/3^(l —then one can 

find that the maximum QFI is Fq^ = 4, which gives the 
minimum value of variance of the estimator S9 = 1/2. 

In the thermodynamics limit, there is no excitation of 
the system in the normal phase, and the atoms part of the 
ground state is the pure Dicke state \j, —j) = Ilfci I 
with al\ Dk = — I i)k- Then the QFI of the reduced 
two-atom state in Eq. (iini) is simplified by FQ{p,a^) = 
4((5cr®)^. For the ground state |j, —j), the variance of 
(Sa^)^ equals to 0 and hence one obtain the minimum 
value of the QFI Fq™ = 0. The minimum value leads to 
the maximum variance of the estimator 50 —>■ oo in the 
normal phase. It also can be easily obtained = 0 

from Eg.f 116 )) with /3 = 0. In the superradiant phase for a 
large coupling strength A Ac with /3 —)> 1/2, the QFI of 
the reduced two-atom state approaches to the maximum 
value Fq^oo —!► 4. In the strong coupling limit, the ground 
state (El) in the Ja;-representation can be described by 

N N 

IG) = ;^[(nje,,)fe)| 0 )A _^/2 + (nj5'^)fe)|0)A^,/2]-where 

CTf|ea;)fe = \ex)k and (Jl\gx)k = -\gx)k- The vacuum co¬ 
herent states | 0 )a±n /2 satisfy A±iv/ 2 | 0 )A±w /2 = 0- It 
yields the variance of (5crf)^ = 1 and the QFI Fq^oo 
equals to 4. It indicates that the variance of the esti¬ 
mator 59 is minimum in the ground state in the superra¬ 
diant phase, which can be distinguished from that in the 
normal phase. 

Fig.([I])(a) displays the QFI of the reduced two-atom 
state in the thermodynamics limit for detunings D = 
A/w = 1. In the normal phase the QFI equals to zero. 
As the the coupling strength enters into the super-radiant 
phase A > Ac, the QFI shows monotonous increasing be¬ 
haviors, demonstrating the existence of the QPT at the 
critical point Ac = VAw/2. As A approaches A —>■ oo 
limit, the QFI Fq^oo tends to the maximum value 4. 


IV. FINITE-SIZE SCALING BEHAVIOR OF 
THE QFI 



FIG. 1: (Color online) QFI Fq^n of the reduced two-atom 
state in the Dicke model as a function of the coupling constant 
A with different sizes N = 20, 256 and N ^ oo for D = 1. (b) 
Scaling of Fq^n as a function of A on a log-log scale at the 
critical point Ac for D — 0.5 and 1, and the solid line scales 

as A-“-«5 ±o-oi. 


gularity coming from (A — Ac)^, one has f{x) ^ x 
which leads to Fq‘//(Ac) ^ , 

In finite-size atom system, the QFI of the reduced two- 
atom state is calculated by the expected value of the col¬ 
lective atomic operators in Eq. (ED in the ground state |G) 
using the bosonic coherent state technique. Eig.([TD(a) 
shows the behavior of Fq^n for different system sizes 
N = 20 and 256. In the normal phase, it tends to zero 
as the number of atoms N increase. And in the super¬ 
radiant phase, A > Ac, the QEI increases from zero to 
the maximum value 4, which are consistent with those 
in the thermodynamics limit. Eor the large system size 
N = 256, there is nearly no deviation of the behavior 
of the QEI from that in the thermodynamic limit. It 
facilitates the calculation of the scaling behavior at the 
critical point by the bosonic coherent state approach. 

We plot the QFI Fq^n as a function of A on a log-log 
scale at the critical point Ac for different detunings D = 
0.5 and 1, as shown in Fig.([TD(b). It is very interesting 
to observe a power law scaling Fq(Ac) oc N'' . Due to the 
advantage of the bosonic coherent state technique [13, 
we are able to study the atom number up to A = 4000 
atoms. The asymptotic slop in the log-log scale for the 
finite size systems gives a exponent v = —0.65 ± 0.01. 

Moreover, the QFI Fq^m in Eq.( [inD can be given ex¬ 
plicitly in terms of the collective atomic operators as 


To study the universality of the superradiant phase 
transition, the exact finite-size scaling law is quite im¬ 
portant. We illustrate the scaling behavior of the QFI 
of the reduced two-atom state. The finite-size scaling 
ansatz for the singular functions Eg*// of the QFI Fq^m 
in the vicinity of the critical point is Cl 

^ ( 17 ) 

where / is a function depending on the scaling variable 
A(A - Ac) 3/2 and n are exponents. To cure the sin- 


32 ( 4 )^ 2N^-8{Jl) 

A(A-l)(A2-2A-b4(J|)) A(A-l) ' 

( 18 ) 

Since the finite-size scaling exponents of the collective 
spin operators in the Dicke model have been derived 
by Vidal and Dusuel 0 , such as (J/)/A^ ^ A , 
(J2)/A2 - A-4/3 and (J2)/A2 - A-^/s. Thus the 
finite-size scaling exponent of the QFI in Eq. (EHD can 
be directly checked. Erom those, it is easily to ob¬ 
tain the finite-size scaling behavior of QFI exactly as 
Fq,n . The result is in consistent with numer¬ 

ical accurate exponent. To the best of our knowledge. 
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such a finite size scaling for the QFI itself has never been 
reported in Dicke model. 

In general quantum metrology process, the QFI de¬ 
pends on the prepared quantum state and the choice 
of the unitary operator. The QFI of the reduced N- 
atom state pA, which is the reduced atomic density ma¬ 
trix by tracing over the field degree of the freedom, has 
been studied in small-size systems TV = 20 in the Dicke 
model [ 2 ^. It is very difficult to predict the finite-size 
scaling exponent due to the too small system sizes investi¬ 
gated based on the basis of the Fock states. Here we focus 
on the exponent of the QFI of the reduced atomic sub¬ 
system for a large number of atoms based on the ground 
state in Eq. (El). The QFI Fa of the reduced atomic state 
Pa{9) = can be given by Eq.(|31), where the 

phase-shift generator U is replaced by Jz and \ipi) are the 
corresponding eigenvectors of pA with nonzero eigenval¬ 
ues {pi}. In the thermodynamics limit, the analytical 
results of the scaled QFI for the atomic subsystem is 


Fa,o./N 


2pA 

e+ +e- + - uj‘^)/{e+ + e_) ’ 


(19) 


where the excitation energies is given by = ^{uj^ + 
A^Z/r^) ± _ A^Z/r^)^ -I- WX'^ujAp. As the cou¬ 

pling strength A approaches to the critical point as A —>■ 
A c, the exc itation energy of tends towards a value of 
\/w^ -f A^ and e_ vanishes, e_ —>■ 0. The critical expo¬ 
nents of the QPT can be manifested in the behavior of 
the scaled QFI 


Fa,c<d 

N 


(A —>■ Ac) 




■ A2 


32w2A? 


A2(16A4-pa;4) 


|Ac-Ap/^ 


( 20 ) 

Fig. (©(a) displays the maximum value 'JuP- F A? jA of 
the scaled QFI Fa^ooIN at the critical point Ac. As ad¬ 
dressed in Ref. [ 2 ^ , the scaled QFI is larger than 1 in the 
normal phase, and then decreases to zero in the superra¬ 
diant phase. 

It is interesting to observe the finite-size scaling be¬ 
havior of the scaled QFI Fa/N in the Dicke model. We 
calculate Fa/N by the ground state [G) in Eq.(El) using 
the bosonic coherent-state technique. Fig.([5])(a) displays 
the scaled QFI Fa/N for finite atom ensemble N = 20 
and N = 256. The scaling behavior of {Fa, oc, — Fa)/N 
as a function of N at the critical point Ac is shown in 
Fig.( [11(b) for different values oi D = 0.5 and 1 on a 
log-log scale. A power-law behavior exists at large N. 
One can see that the finite-size exponents extracted from 
all curves tend to a converging value —0.33 ± O.OI in 
Fig.(|l])(b). The precise estimate of the exponent for the 
QFI is very significant to help clarify the universality of 
the QPT. 




FIG. 2: (Color online) (a)The scaled QFI Fa/N of the re¬ 
duced A-atom state as a function of A with different sizes 
N = 20, 256 and A —>■ 00 for D — 1. (b)Scaling of the scaled 
QFI {Fa,00 — Fa)/N as a function of A on a log-log scale at 
the critical point Ac for D = 0.5 and 1, and the solid line 
scales as iV-o-33±o.oi^ 


QPT in the Dicke model from the quantum estimation 
perspective. The QFI are obtained analytically in the 
thermodynamics limit. The behavior of the QFI of the 
reduced two-atom state shows that the ground state un¬ 
dergoes a superradiant phase transition at the critical 
regime in the infinite size system. And Finite-size scal¬ 
ing exponents of the QFI are calculated up to a large 
atom number N = 4000 by the bosonic coherent tech¬ 
nique. Power law scaling behavior at the critical point 
is observed. Moreover, large-A scaling behavior of the 
scaled QFI of the reduced A-atom state is also calculated, 
giving the accurate exponent. Such a scaling behavior of 
QFI has not been reported in the critical systems of Dicke 
model, as far as we know. These salient features might 
be used for quantum metrology and quantum estimation 
in some experimentally realized systems to the quantum 
information science and the quantum computing. 
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Appendix 


V. CONCLUSION 


For the reduced two-atom state p in Eq. (0), the cor¬ 
responding eigenvalues are given by 


In summary, we have proposed the QFI of the reduced 
two-atom state to capture the quantum criticality of the 


Pi = 2 w,P2 = 0,p± = -{v+ -f u_ ± v^). 


( 21 ) 
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where 7 = — f_)^ + 4|up. The corresponding eigen¬ 

states to Pi and p± are 


l<^±) 




( 22 ) 


(23) 


with = 2uj ± (t!+ — For the unitary oper¬ 

ation S = exp(—), the QFI of the reduced two-atom 


state can be evaluated as 


Fq = ^P±{5al)\ +‘ipi{5al)'{- 


16 p+p- 
p+ +p- 


-E 

i=ih 


ISpiPi 

Pi +pi 




(24) 


where the variance of operator erf is (derf)! = 
((/)i|(crf)^|((>i) — \{(j)i\a^^\(j)i)\'^. Substituting the values of 
P±,i and \4>±,i) into the above equation, the QFI can be 
given analytically. 
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